Effects of symmetry energy in the reaction ''°Ca+^^''Sn@140 MeV/nucleon 
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The density dependent symmetry energy is a hot topic in nuclear physics. Many laboratories in 
the word are planning to do related experiments to probe the symmetry energy. Based on the semi- 
classical Boltzmann-Uehling-Uhlenbeck (BUU) transport model, we studied the efTects of nuclear 
symmetry energy in the central reaction '"'Ca+^^^Sn at 140 MeV/nucleon in the laboratory system. 
It is found that the rapidity distribution of free nucleon's neutron to proton ratio is sensitive to the 
symmetry energy, especially at large rapidities. The free n/p ratios at small or large rapidities may 
reflect high or low density behavior of nuclear symmetry energy. To probe the density dependence 
of nuclear symmetry energy, it is better to give the kinetic distribution and the rapidity distribution 
of emitted nucleons at the same time. 

PACS numbers: 25.70.-z, 21.65.Ef 



I. INTRODUCTION 

After about three decades of intensive efforts in both 
nuclear experiments and theories, the equation of state 
(EOS) of isospin symmetric nuclear matter is now rela- 
tively well determined mainly by studying collective flows 
in heavy-ion collisions and nuclear giant monopole reso- 
nances [l|, 0. The major remaining uncertainty about 
the EOS of symmetric nuclear matter is due to our poor 
knowledge about the density dependence of the nuclear 
symmetry energy [H, 0-13|- Therefore, the new challenge 
is to determine the EOS of isospin asymmetric nuclear 
matter, especially the density dependence of the nuclear 
symmetry energy. Besides the great interests in nuclear 
physics, the EOS of asymmetric nuclear matter is also 
crucial in many astrophysical processes, especially in con- 
nection with the structure of neutron stars and the dy- 
namical evolution of proton- neutron stars [ll|. Consid- 
erable progress has been made recently in determining 
the density dependence of the nuclear symmetry energy 
around the normal nuclear matter density. However, 
much more work is still needed to probe the high-density 
behavior of the nuclear symmetry energy. Currently, 
to pin down the symmetry energy, the National Su- 
perconducting Cyclotron Laboratory (NSCL) at Michi- 
gan State University, the Gesellschaft fuer Schwerioncn- 
forschung (GSI) at Darmstadt, the Rikagaku Kenkyusho 
(RIKEN, The Institute of Physical and Chemical Re- 
search) of Japan, and the Cooler Storage Ring (GSR) 
in Lanzhou, are planning to do related experiments. Ef- 
fects of the symmetry energy on the pre-equilibrium neu- 
tron/proton ratio in heavy ion collisions have been ana- 
lyzed with many authors. For example, calculations by 
Zhang et aL^Sj] within the Improved Quantum Molecular 
Dynamics(ImQMD) model using two different density- 
dependent symmetry-energy functions, the results of neu- 
tron to proton ratios are sensitive to the density depen- 
dence of the symmetry energy. Within the transport 
model Lanzhou Quantum Molecular Dynamics(LQMD), 



Feng investigated the single and double neutron to 
proton ratio of free nucleons with different collision cen- 
tralities and effective mass splitting. Recently, Ma et 
al.[Toj studied the high-density behavior of the symme- 
try energy by using single and double ratios of neutrons 
to protons within Isospin Quantum Molecular Dynam- 
ics(IQMD) model, they confirmed that it is possible to 
study the high-density behavior of symmetry energy by 
using the neutron-to-proton ratio from free nucleons. In 
the present work we give our results of the effects of sym- 
metry energy on free neutron to proton ratio as a func- 
tion of rapidity in the '^°Ca-h^^'*Sn at 140 MeV/nucleon. 
We find that nucleon emissions at low-angle and large 
rapidities regions are suitable to be used to probe the ef- 
fect of symmetry energy. The free n/p ratios at large or 
small rapidities may reflect low or high density behavior 
of symmetry energy. 



II. THE IBUU04 TRANSPORT MODEL 

Our present study is based on the transport model 
IBUU04. In this semi-classical model besides nucleons, 
A and N* resonances as well as pions and their isospin- 
dependent dynamics are included, the experimental free- 
space nucleon-nucleon (NN) scattering cross sections and 
the in-medium NN cross sections can be used option- 
ally. By using the relativistic mean field theory, the 
initial neutron and proton density distributions of the 
projectile and target are obtained. The isospin depen- 
dent phase-space distribution functions of the particles 
involved are solved by using the test-particle method nu- 
merically. The isospin-dependence of Pauli blockings for 
fermions is also considered. More details can be found in 
Refs. [T2I4T7I. In the present work, we use the isospin- 
dependent in-medium NN elastic cross sections from the 
scaling model according to nucleon effective masses jl3j . 
In the IBUU04 transport model, the most important in- 
put is the momentum- and isospin-dependent single nu- 
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cleon potential, we use a single nucleon potential de- 
rived within the Hartree-Fock approach using a mod- 
ified Gogny effective interaction [l^, the momentum- 
dependent single nucleon potential (MDI) adopted here 
is: 



U{p,6,p,T) = Au{x)^ + Mx)^ 
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In the above equation, S = {pn — Pp)/{Pn + Pp) is the 
isospin asymmetry parameter, p = Pn + Pp is the baryon 
density and p„ , pp are the neutron and proton densi- 
ties, respectively, t = l/2(— 1/2) for neutron (pro- 
ton) and T ^ t' , a = 4/3, /r(r, p) is the phase-space 
distribution function at coordinate r and momentum 
p. The parameters Au{x), Ai{x), B,Cr,T, Cr,T' and A 
were set by reproducing the momentum-dependent po- 
tential U{p, S, p, r) predicted by the Gogny Hartree-Fock 
and/or the Brueckner-Hartree-Fock calculations. The 
momentum-dependence of the symmetry potential stems 
from the different interaction strength parameters Cr,T' 
and C'r,T for a nucleon of isospin t interacting, respec- 
tively, with unlike and like nucleons in the background 
fields, more specifically, Cuniike = —103.4 MeV while 
Ciike — —11.7 MeV. The parameters yl„(x)and Ai{x) de- 
pend on the x parameter according to Au{x) — —95.98 — 
x-§^ and A/(x) = -120. 57-1- a; The saturation prop- 
erties of symmetric nuclear matter and the symmetry en- 
ergy of about 32 MeV at normal nuclear matter density 
Po = 0.16 fm~^. The incompressibility of symmetric nu- 
clear matter at normal density is set to be 211 MeV. Ac- 
cording to essentially all microscopic model calculations, 
the EOS for isospin asymmetric nuclear matter can be 
expressed as 



E{p, S) = E{p, 0) + E,y„,{p)S^ + 0(6'^), 



(2) 



where E(p, 0) and i?sym(p) are the energy per nucleon of 
symmetric nuclear matter and nuclear symmetry energy, 
respectively. For a given value x, with the single parti- 
cle potential U{p,S,p,T), one can readily calculate the 
symmetry energy i?sym(p) as a function of density. 



III. RESULTS AND DISCUSSIONS 

Fig. [1] shows the density dependence of nuclear sym- 
metry energy with parameter x — 1,0, respectively. As 
discussed in the previous part, the single particle po- 
tential used has an x parameter, different specific x 
parameter denotes different density dependent symme- 
try energy. For the central reaction '*''Ca-|-^^'*Sn at 140 
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FIG. 1: (Color online) Density dependence of nuclear sym- 
metry energy with parameters a; = 1,0, respectively. 
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FIG. 2: (Color online) Maximal baryon density reached in the 
central reaction ""Ca-l-^^^Sn at 140 MeV/nucleon. 



MeV/nucleon, the maximal density reached is about 1.5 
times saturation density as shown in Fig. [5] We there- 
fore only show the low density's symmetry energy as a 
function of density. From Fig. [1] we can also see that the 
low density behaviors of nuclear symmetry energy sepa- 
rate from each other with different x parameters. At the 
saturation point there is a cross and then they separate 
from each other again. At lower densities, the value of 
symmetry energy of a; = is lower than that of a; = 1. 
But at high densities, the value of symmetry energy of 
x = is higher than that of x = 1. From Fig. [51 we can 
see that at the reaction time i = 50 fm/c, the reaction 
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almost ends. 
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effects of symmetry energy are always very small. The 
maximal longitudinal momentum at positive momentum 
part reached is about 0.9 GeV/c, but the maximal longi- 
tudinal momentum at negative momentum part reached 
is only about 0.3 GeV/c. From this plot we can clearly 
see that most emissions are at very small longitudinal 
momentum (pz), denoting that more nucleons may emit- 
ted perpendicularly to the reaction plane. Experimen- 
tally, one always needs to know the angle distribution 
of the emitted number of probed nucleons in the reac- 
tion. This is because one usually can not probe the emit- 
ted particles in At: directions. For this purpose, we plot 
Fig. |4l the angle distribution of relative emitted number 
{N cos 9 /Ntotai) of the free nucleons in laboratory system 
in the central reaction '"'Ca-|-^^'*Sn at 140 MeV/nucleon. 
From this plot, we can clearly see that whether for x = 1 
or X = nucleons are inclined to low-angle emission. It 
is thus suitable to probe the emitted nucleons in the low- 
angle regions. 



FIG. 3: (Color online) Longitudinal momentum (p^) distri- 
bution of emitted free nucleons in laboratory system in the 
central reaction "^^Ca+^^^Sn at 140 MeV/nucleon. 
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FIG. 5: (Color online) Neutron to proton ratio (n/p) of the 
emitted free nucleons as a function of reduced rapidity with 
X — 1 and in laboratory system in the central reaction 
''°Ca+^2^Sn at 140 MeV/nucleon. 



FIG. 4: (Color online) Angle distribution of relative emitted 
number {N cos 8 / Ntotai) of the free nucleons in laboratory sys- 
tem in the central reaction '"'Ca-|-^^'*Sn at 140 MeV/nucleon. 

Fig. |3] shows the longitudinal momentum (^2) distri- 
bution of emitted free nucleons in laboratory system in 
the central reaction '"^Ca-l-^^^Sn at 140 MeV/nucleon. In 
our work, free nucleons are identified as those having lo- 
cal baryon densities less than p — ipo- We can first see 
that more nucleons are emitted at positive momentum 
part than at the negative part. Also one can see that 
longitudinal momentum (pz) distribution of emitted free 
nucleons is insensitive to the symmetry energy since the 



To see the effects of symmetry energy on emitted free 
neutron to proton ratio, we provide Fig. [SJ the neutron 
to proton ratio (n/p) of the emitted free nucleons as a 
function of reduced rapidity with a; = I and in labo- 
ratory system in the central reaction ^°Ca-|-^^^Sn at 140 
MeV/nucleon. From this plot we can see that at large 
rapidities, effects of symmetry energy are clearly shown 
on the neutron to proton ratio. The soft symmetry en- 
ergy (x = 1) corresponds large n/p ratio than the stiff 
case {x = 0), the n/p ratio thus reflects the low-density 
behavior of nuclear symmetry energy. Around low rapid- 
ity y/ybeam= 0.2, the Stiff Symmetry energy corresponds 
large n/p ratio, reflecting the high-density behavior of 
nuclear symmetry energy. This indicates that nucleons 
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FIG. 6: (Color online) Neutron to proton ratio (n/p) of the 
emitted free nucleons as a function of kinetic energy with x= 1 
and in laboratory system in the central reaction '"'Ca+^'^'^Sn 
at 140 MeV/nucleon. 



around low rapidity y/ybeam^ 0.2 are from high-density 
region of compressed nuclear matter, and this part of 
emitted nucleons has a large number (also as shown in 
Fig. 131). Thus the free n/p ratios at large or small rapidi- 
ties may reflect low or high density behavior of symmetry 
energy. And because the emitted nucleons are mainly at 
low-angle regions (shown in Fig. 2]), we deduce that the 
better way to probe the effects of symmetry energy on 
free neutron to proton ratio n/p is probing the emitted 
nucleons at large rapidities and low-angle regions in lab- 
oratory system in the central reaction '^^Ca-f ^^^Sn at 140 
MeV/nucleon. 

Fig. [6] shows the neutron to proton ratio {n/p) of the 
emitted free nucleons as a function of kinetic energy with 
X = 1 and in laboratory system in the central reac- 
tion 4°Ca+i24gj^ MeV/nucleon. We can see that 
the stiff symmetry energy [x — 0) corresponds large n/p 
at higher kinetic energies. At lower kinetic energies the 
soft symmetry energy {x = 1) corresponds larger n/p. 
From Fig. [6] and Fig. [5l we can again deduce that nucle- 
ons at small rapidities mainly come from the compressed 
nuclear matter, which have large kinetic energies, their 
n/p ratios reflect high density behavior of nuclear sym- 
metry energy. Nucleons at large rapidities mainly come 
from the low-density nuclear matter, which have no very 
large kinetic energies, their n/p ratios reflect low den- 
sity behavior of nuclear symmetry energy. From above 



discussions, in the central reaction ^''Ca-|--'^^''^Sn at 140 
MeV/nucleon, to probe the high density behavior of nu- 
clear symmetry energy, we can probe nucleons at very 
high kinetic energies and at mid-rapidities. To probe 
the low density behavior of nuclear symmetry energy, we 
can probe nucleons at not very high kinetic energies and 
at large rapidities. Therefore, to probe the symmetry 
energy, it is better to give the kinetic distribution and 
the rapidity distribution of emitted nucleons at the same 
time. 



IV. SUMMARY 

Based on the semiclassical Boltzmann-Uehling- 
Uhlenbeck (BUU) transport model, we studied the 
effects of symmetry energy on the free neutron to 
proton ratio in the central reaction ^°Ca-|-^^''Sn at 140 
MeV/nucleon in the laboratory system. It is found 
that at large rapidities free nucleon's neutron to proton 
ratio is quite sensitive to the symmetry energy. And 
we also find that free nucleons are mainly emitted at 
low-angle regions and at zero longitudinal momentum 
(pz) in the laboratory system. The maximal longitu- 
dinal momentum at positive momentum part reached 
is about 0.9 GeV/c, but the maximal longitudinal 
momentum at negative momentum part reached is 
only about 0.3 GeV/c. And more nucleons are emitted 
at positive momentum part (mainly in the range of 
Q < Pz < O.SGeV'/c) than at the negative part. The free 
n/p ratios at small or large rapidities may reflect high or 
low density behavior of nuclear symmetry energy. These 
information is useful for related experiments relevant to 
symmetry energy studies. 
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